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Within an effective hadronic model including electromagnetic interactions via a Uem(l) gauge, we 
reinvestigate photon Bremsstrahlung from a hot hadronic gas as expected to be formed in relativistic 
heavy-ion collisions at SPS energies. We calculate photon emission from the reactions tttt — > nn^ 
and -kK — > ■nK'y by an explicit (numerical) evaluation of the multi-dimensional phase space integral. 
This, in particular, allows to avoid the commonly employed soft photon approximation (SPA), as 
well as to incorporate final-state thermal enhancement factors. Both improvements are shown to 
result in an appreciable increase of the photon production rate over previous hadronic calculations. 
Upon convolution over a thermal fireball we find an improvement in the description of recent low 
transverse- momentum WA98 data at SPS. The influence of both Landau- Pomeranchuk-Migdal and 
f~^ ' in-medium effects on "ct" and p-meson exchanges are briefly discussed. 

o' 

f^^ ' PACS numbers: relativistic heavy-ion collisions, photon spectra, Bremsstrahlung 
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5j ■ I. INTRODUCTION 



A hot deconfined state of matter called quark-gluon plasma (QGP) is predicted to have existed during the first few 

microseconds after the big bang. Ultrarelativistic heavy-ion collisions (URHICs) provide the only way to recreate such 

CO matter for a short moment in the laboratory. Electromagnetic radiation is expected to be a unique probe in the study 

K^ of URHICs, as the (real and virtual) photons, once produced, decouple from the strongly interacting hot matter. 

^~; ■ Since photons are emitted throughout the entire evolution of the fireball, including QGP and hadronic phases, their 

J^ \ spectra carry valuable information on the various environments of their creation. 

— ^ I Studies of photon emission from hot and dense matter in heavy-ion collisions have a long history both theoretically 

^) I and experimentally, cf., e.g., Refs. [ll, 0jB0, Q for recent reviews. While scattering in deconfined matter is expected 
\^ , to mostly manifest itself in the photon spectra at transverse momenta above qt — 2 GeV, the hadronic stages are likely 
^^ ' to dominate the emission at low qt [y, 0, H, M, llO) llll) E) U3\ ■ The latter regime has recently received renewed interest 
due to measurements of the WA98 collaboration in central Pb-Pb collisions at the CERN Super Proton Synchrotron 
(SPS) [ij,[i5|. Using Hanbury Brown- Twiss (HBT) interferometry methods, a direct photon signal could be extracted 
in the range of gt — 0.1 — 0.3 GeV 'l5l|. These data exhibit a large excess over theoretical predictions [lil, [i3| that have 
previously shown agreement with data for the same system (and by the same experiment) 14] at higher momenta, 
qt > 1.5 GeV (where the more conventional subtraction method enabled the extraction of a direct photon excess). In 
subsequent work [i^|lajll3|! ^^^ ^'^^^ of photon Bremsstrahlung from elastic tt-tt interactions, which was not included 
in the previous calculations, has been re-assessed. It was found that the Bremsstrahlung contribution notably increases 
the spectral yield at low qt by up to 50% over the baseline predictions, but still significantly falls short of the WA98 
data [13. 

Cu I The main objective of the present work is to revisit thermal Bremsstrahlung from elastic meson-meson interactions. 

Based on an effective hadronic Lagrangian, with electromagnetic interactions implemented via Ucm(l) gauging, we im- 
prove the accuracy of earlier analyses by going beyond the commonly employed soft photon approximation (SPA) [18| , 
and extend the calculations to the strangeness sector (jr-K scattering). In the context of soft dileptons, corrections to 
the SPA have been found to be sizable and positive [l9|. Our approach furthermore enables the explicit inclusion of 
thermal final-state enhancement factors, as well as meson-chemical potentials which arise in the hadronic evolution 
of URHICs after hadrochemical freezeout. As an indication of the uncertainties in applications to heavy- ion data, we 
briefly discuss other meson- meson Bremsstrahlung sources 112|, as well as schematic estimates of medium effects in 
the scalar tt-tt channel (dropping a mass) and of the Landau-Pomeranchuk-Migdal (LPM) effect [20, UM ■ The tt-tt 
and TT-K Bremsstrahlung rates are then combined with earlier hadronic emission calculations [ll| and convoluted over 
a thermal fireball model for central Pb-Pb collisions at SPS in order to put the improvements into context with the 
recent low- momentum WA98 data flSj . 

Our article is organized as follows. In Sec. |TT] we introduce the chiral effective Lagrangian that will be used to 
evaluate amplitudes for elastic tt-tt and n-K interactions, with parameters fixed to reproduce the corresponding vacuum 
scattering data in Born approximation. In Sec. IIIIl we discuss our treatment of the basic kinetic theory expression to 
calculate Bremsstrahlung from tt-tt and n-K scattering in a thermal bath. We first construct the pertinent scattering 
amplitudes by gauging the hadronic Lagrangian (Sec. IIII A[) . with special care to maintain electromagnetic gauge 
invariance by implementing appropriate contact interactions (4-point vertices). We then give a detailed account of 
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the various features that govern the emission rate (Sec. IIIIB[) . with emphasis on the improvements over previous 
evaluations. In Sec. IIVI we illustrate further uncertainties when applying low-energy rates to heavy-ion collisions in 
terms of schematic estimates of the LPM effect (Sec. IIV A|) and medium modifications of the exchanged p and a 
mesons (Sec. IIV Bj) . In Sec.|Vl the soft photon production rate is applied to calculate photon spectra in central Pb-Pb 
collisions at SPS, which arc supplemented with earlier calculations and compared to experimental data from WA98. 
In Sec. IVI[ we estimate the coherent photon emission from initial scattering of protons within both projectile and 
target nuclei. Sec. IVIIJ contains a summary and discussion. 



II. ELASTIC MESON-MESON SCATTERING 

In this work we focus on photon Bremsstrahlung induced by elastic interactions of the lightest constituents in a 
hadronic gas, i.e., pion-pion (tt-tt) and pion-kaon (tt-K) scattering (the latter will turn out to produce a factor of 4-5 
less low-energy photons than the former; other contributions will be briefly discussed at the end of Sec. IIIIBp . We 
first recall the underlying model 13, '22\ that will allow us reproduce vacuum scattering data and thus form the basis 
to evaluate the pertinent Bremsstrahlung processes, tt^tt'' -^ tt^tt'^j and ■n°'K^ -^ n^K'^j. Starting point is a chiral 
Lagrangian for the pseudoscalar fields (tt and K) [23| where the low-lying vector mesons, p and K* , are introduced 
via a covariant derivative leading to 



-igrrKK'iKfK*^ 



Sp-rniP 



n X df,n) + igpKKiKrdf^K 
d^.KrK*'^ ■ 7?) + h.c. . 



d^KfK) ■ ^ 



(1) 



Since we will restrict ourselves to the Born approximation, an additional (chirally invariant) cr- tt-tt interaction vertex 
with a low-mass "cr" -meson ]2A] has been added to adequately account for S'-wave tt-tt scattering near threshold. The 
corresponding Feynman diagrams, comprising s, t and u channel exchanges, are depicted in Fig. [T] The pertinent 
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FIG. 1: Diagrams for tt^tt*' -^ ■k^-k'^ scattering. 




amplitudes for elastic tt^-tt scattering in t- and s-channel are given by 

,. , . -4:gl^^F^{t){pa-pi){pb-p2) , ~gl^^F^it){Pa-Pb+Pa-P2+Pb-Pl+Pl-P2) 

Mt{Pa,Pb,Pi,P2) = : :r—. t; + 



t-m'i+ iniaVo 



t -mf, + irUpTp 



Ms{Pa,Pb,Pl,P2) 



'Agl^^{pa ■Pb){pi ■P2) g^TrAPa ' Pi - Pa ' P2 + Pb ' P2 ~ Pb ' Pi) 



inicrTa- 



inioTo 



(2) 



To account for the finite size of the hadronic vertices, we introduce a momentum-transfer damping monopole form- 
factor, F(q^), for the t and u channels [9i. i22i]. 



Fiq' 



^l-q^ 



(3) 



where ma is the mass of the exchanged meson and q the four-momentum transfer (for simplicity, we neglect formfactors 
for s-channel polegraphs). The differential cross section for tt+tt^ — > tt+tt^ then follows as 



da _ \Mt+Ms\^ 
dt 64TTspf„ 



(4) 



with Pin = ^/{s — Am'^)/2 the initial pion three-momentum in the center-of-mass (cm) frame. Fixing the parameters 
in the non-strange sector as in Refs. [a, [24I, m^ = 0.525 GeV, F^ = 0.1 GeV, g^j^^ma- = 1.85, rup = 0.775 GeV, 



Tp — 0.155 GeV and ^ptttt = 6.15, provides a fair fit to the total tttt cross section [25| up to a cm energy of ^1 GeV, 
cf. left panel of Fig. [51 We neglect the contribution from the /2(1270) tensor meson as its thermal density is a factor 
of ^--^25 smaller than p-density. In the right panel, we display the decomposition into a- and p-meson exchanges which 
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FIG. 2: Left panel: data on the elastic ■k'^tv — > tv^tt cross section [25| compared to a fit using a meson-exchange model in 
Born approximation. Right panel: decomposition of the cross section into contributions from a- and p-exchange. 

identifies the s-channel p-polegraph as the prevalent contribution to elastic tt^tt^ scattering, while the one from 
cr-exchange is small. 

Next, we turn to tt-K scattering, Bremsstrahlung off which could be significant due to both the kaon's relatively 
large abundance in the hadronic fireball and its resonance cross section via the K* resonance. The relevant Feynman 
diagrams follow from Fig.[T]by replacing one of each in- and outgoing pion by a kaon, and the p by K* (no cr-exchange), 
resulting in sixteen processes with charged particles in the initial and/or final state. The corresponding amplitudes 
for elastic Tr~-K'^ scattering are given by 



-gp7r-n-gpKKF'^it)iPa ' Pb + Pa ' P2 + Pb ' Pi + Pi ' P2) 



MtiPa,Pb,Pl,P2) = 
Ms{Pa,Pb,Pl,P2) = 



As before, the hadronic t- and w-channel vertices are augmented by a formfactor. 
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m\. 





Al-q^ ' 



(5) 



(6) 



which is of the same form as above but with a slightly different parameter dependence which facilitates a better fit 
to the free scattering data. Applying Eq. ^ and fixing A = 1.5 GeV, Tk* — 0.051 GeV, gpKK = g-nKK* = 3.55, and 
rriK' — 0.89 GeV, the empirical cross section for elastic -k^ K^ -^ tt^ K^ scattering [2y| is reasonably well described, 
cf. Fig. H 



III. THERMAL PHOTON EMISSION RATE 



In kinetic theory, and to leading order in the e.m. coupling constant a^m = 1/137, the thermal emission rate (per 
unit four- volume) of photons of energy qo = q for a Bremsstrahlung process of type a + b —> 1 + 2 + 7 can be cast into 
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FIG. 3: Total elastic tv K^ ^ tt K cross section as calculated in a meson-exchange model in Born approximation compared 
to experimental data from Ref. [2q |. 



the form 



qo 



ClR-y 



= N 



(fipa cPpb (Ppi d^P2 

2Ea{2'Kf 2Eb(2T:f 2Ei{2nf 2E2{2-Kf 



X {2T:Y5'^{pa+Pb - Pi - P2 - g)|A^7l' 



. J{E,)f{Eb)[l±f{E^)][l±f[E2)] 
2(27r)3 



(7) 



where M.^ is the corresponding scattering amplitude, f{Ei) = l/(e(-^'~^*)/^ ± 1) = /^ are Fermi-Dirac or Bose- 
Einstein distribution functions {fif. chemical potential of hadron i), and the overall degeneracy factor TV depends on 
the specific process. As elaborated in Appendix El we rewrite the emission rate without further approximation as 



qo- 



dR^ 



N 



d^q 16(27r)ii 

where A is given by 



IpJdEadq^a I \pb\dEbSme,,d9bd<l)b f \pi\dEism9id9id<Pifafb[l + /i][l + /2 



\M, 



A 



A^\^'{cosei)\ 

and cos6'!j is the root of the function (p{cos9ii), 

ip{cos9a) = (2|pQ||pf,|sin6'bCOS0aCOS0b + 2|pQ||pf,|sin6'bSin<^aSin(/)b - 2|pa||pi|sin6'icos0acos0i 
-2|pa| |pi |sin6'isin(?!)asin(/)i)^l - cos^^a 

+ (2|pa||p6|cOs6'b - 2|pa||pi|cOS0i - 2|pa| |q| )cOs6'a 
+ |PaP + iPb - Pi - q)' - (Ea +Eb-Ei- qof 



(8) 



(9) 



mi =0 



(10) 



The main point of the present work is that the numerical evaluation of Eq. ([5]) goes beyond previous treatments 0, [24I 
by including the final state enhancement factors, (1 + /i) x (1 + /2)^, and by evaluating Ai^ beyond the soft photon 
approximation (SPA). The latter point is detailed in the following Section, where a summation over all possible 
processes for thermal photon Brcmsstrahlung from vr-vr and t:-K scattering is performed. 



^ The relevance of the final-state enhancement factors in hadronic phase of URHICs is further augmented by the build-up of substantial 
pion- and kaon-chemical potentials after chemical freezeout [2711 . 



A. Amplitudes for Bremsstrahlung 

The photon coupHng to pseudoscalar and vector mesons can be implemented via a Uem(l) gauge [7|], rendering the 
interaction Lagrangian with electromagnetism as 

£ = eA^'{^f,^l x 7r)3 + 2eg^^^crA''(a^7r x 7r)3 - egp^^A^,[K x (7? x p^)]3 

+e{A>^{dp(r X p,)3 + [{d^A'^p, - A'^d^p,) X p^]3 + [p^ X {A'^d^p, - d^A^p^M 
+ie[A>'{K*''QdpK: - dpK*''QK:) + K*>'Q{d^A'' Kl - A'^d^Kl) + {A'^B^Kl - dA''K:)QK*^] 
+ieA>'{KQd^K - df,KQK) - eg^KK-A^'[k*{2fQ - Qf)K + K{2Qf - tQ)K*] ■ 7? 
+egpKKA^K{fQ + Qf)K-p^,, (11) 

where A^^ is the electromagnetic field, Q=diag(l, 0) the charge operator, and the subscript "3" denotes the third 
component of isospin. 

First, we consider the process for Bremsstrahlung of a photon via tt^tt — > 7r^7r^7 scattering, based on the diagrams 
depicted in Fig. [TJ In charge basis, this amounts to the following seven processes: tt+tt" -^ 7r+7r"(7r*'7r*'), 7r+7r° — > 
TT+vr*', TT^TT*^ -^ tt^tt'', tt'^tt'' — > TT^'^Tr^, 7r+7r+ — !■ 7r+7r+ and tt^tt^ — > tt^tt^, with a photon attached to each external 
or internal charged particle. E.g., the amplitude for the process tt+tt^ — > 7r+7r^7 with a photon attached to a charged 
external pion in all possible ways reads 

M^ ^ eJ^a[-Mt{Pa-q,Pb,Pl,P2)+Ms{Pa-q,Pb,Pl,P2)] 

+ eJl^[Mt{Pa,Pb - q,Pl,P2) + Ms{Pa,Pb - q,Pl,P2)] 
+ ejf [Alt(pa,P6,Pl + q,P2) + Ms{Pa,Pb,Pl + q,P2)] 

+eJ^[Mt{pa,Pb,pi,P2 + q) + Ms{pa,Pb,pi,P2 + q)] , (12) 

where Aig and Ait are the amplitudes quoted in Sec. [ill In the above, 

j^ _ -Qa.bi^Pa.b - q)" J^, _ Qi,2(2pi,2 + g)^ .^„. 

'''"- 2pa.b-q ' ^'^" 2p,.,-q ^^'^' 

are the electromagnetic currents with Qi the charge of the pion in units of the proton charge. Electromagnetic 
gauge invariance furthermore requires the incorporation of contact diagrams with the photon attached to each proper 
vertex. In Appendix iBl we explicitly construct the pertinent amplitudes so that the total satisfies q^ ■ M'^ = 0. For 
the processes with four-point vertices one arrives at 



Mti 



^egl^^F[{pb - P2f]{pa + PiTPb ■ P2 ^egl^^F[{pa ~ pif]pa ■ pi{pb + P2T 
ipb - P2Y - ml + irriaTcr (pa - Pi)'^ - ml + im^T^ 

4e5^^„ {Pa -PbYpi-P2 Aegl^^pa ■ Pb{pi - P2)'' 



{pi +P2)^ - ml + imaTa [pa + PbY -m1 + im^Ta 
l%.F[{pa-prf]{pa+p^Y 2egl^, 



2egl^^F[{pa-pif]{Pa+PiY , 2egl^,F[{pb - p2f]{pb + P2Y 



(Pa - PiY - m-p + impTp {pb - P2Y " "^p + impTp 



'^egL^Pa - PbY 2eg2 (p^ „ p^Y 



{Pa + PbY - mj, + impFp (pi + P2Y - m?p + impFp 



(14) 



The same procedure of deriving the Bremsstrahlung amplitude, leading to Eqs. (|12p . can be applied to tt-K scat- 
tering. The additional contact terms necessary to maintain gauge invariance are given by 

j^p ^ '^e9p7,T^9pKKF[{pa-pif](pa+piY ^egp^^^gpRK F[{pb - P2Y]{Pb + P2Y 
{Pa~ PiY -m^p + impTp [pb - P2Y - m'^ + impTp 

^eglKK-AiPa-PbY - {rnl - mg)(p„ +pfc)^/m|..] 



{Pa + PbY - m\, + imK-Tx' 
iIkk'UPi -P2Y-iml- ml){pi +P2YI'm\.] 
{pi + P2Y - ^K" + imx'TK' 



(15) 



which completes the set of processes for ttK -^ ttK^. 
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FIG. 4: Comparison of various approximations to the thermal photon rate for tttv -^ tvtt^ reactions as a function of photon 
energy. 

B. Soft Photon Emission Rates 



With the ampHtudes as specified in the previous section we now proceed to the results of the numerical integration 
of the thermal photon production rate, Eq. ([8|). In Fig. |4] we compare our new results for Brcnisstrahlung off tt-tt 
scattering to earlier calculations. With the application to URHICs in mind, we choose medium conditions representing 
Pb-Pb collisions at full SPS energy, roughly half way through the hadronic evolution, corresponding to a temperature 
of T = 150 MeV and a pion chemical potential of /^^r — 40 MeV [ll].l28l|. When working in SPA and neglecting final- 
state Bose enhancement factors for the outgoing pions, our results closely coincide with computations by Turbide et 
al. [i^[l3] (dashed vs. solid line in Fig. |4]). When the SPA is relaxed, the rate is found to increase by about 25% at 
low energies (go = 0.1 GeV) and close to 40% at higher energies {qq — 0.4 GeV), see the short-dashed line in Fig. U) 
Finally, when implementing the final-state Bose enhancement, the rate rises another ^30%. Thus, the total increase 
of soft photon emission over previous computations amounts to 60-100%, which is quite appreciable and the main 
result of our paper. 

In Fig. [5] we display the temperature evolution of the Bremsstrahlung rate from tt-tt scattering. The left panel 
includes pion-chemical potentials as estimated in Ref. [28l | to preserve the pion (and kaon) number in the hadronic 
evolution of Pb-Pb collisions after chemical freezeout, parametrized by a linear increase with temperature (T in units 
of GeV), 



^- ^ ^:o^ ^^^ 



AiK = 1.9 (0.175 -T) GeV 



(16) 



(a small nonvanishing value of fi.^ at chemical freezeout has been introduced to reproduce the measured pion-to- 
baryon ratio at full SPS energy; this is due to the somewhat limited number of resonances included in the equation 
of state in Ref. ^] which affects the pion yield via feeddown; it slightly overestimates (underestimates) the direct 
pion density at large (small) temperature, e.g., e^"/-^ = 1.11 at T = 0.18 GeV). When reducing the temperature from 
T = 0.18 GeV to 0.11 GeV, the decrease in rate is close to a factor of 10 at qo = 0.1 GeV but becomes a factor 
of ~15 at qo = 0.4 GeV, reflecting the larger slope at smaller T. This also entails an increasing weight of photon 
emission at lower energies toward later stages in the fireball evolution (note that increase in the hadronic fireball 
volume between chemical {Tc — 0.175 GeV) and thermal freezeout {T ~ 0.110 GeV) is about a factor of ^6, implying 
that the yield of Bremsstrahlung photons is rather sensitive to the late stages of the fireball lifetime) . In the right 
panel of Fig. [3 the pion-chemical potentials are set to zero: the decrease compared to the ^^ ^ results amounts to 
slightly more than a fugacity factor squared, e^^^'"^ (as one would naively expect for tttt annihilation), mostly due to 
the final-state enhancement factors. Finally, a comparison to the hadronic rate calculations of Ref. [11| reveals that 
TT-TT Bremsstrahlung as computed here exceeds previously dominant baryonic contributions by close to a factor of 2 
at go = 0.2 GeV, and is smaller by a factor of 2 at go = 0.4 GeV. 
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FIG. 5: Thermal photon emission rate from Bremsstrahlung via mi -^ mvf reactions as a function of photon energy at different 
temperatures. The left panel includes pion chemical potentials appropriate for P6(158 AGeV)-P6 collisions while in the right 
panel fj.^ = 0. 
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FIG. 6: Thermal photon Bremsstrahlung emission rate from irK — > nKy reactions as a function of photon energy at different 
temperatures. The left panel includes pion- and kaon-chemical potentials according to Eq. p6p . while in the right panel 

fJ-TT = fJ-K = 0. 



Turning to 7 radiation off elastic w-K scattering (Fig. [S]), we find the emission rate to be at the ~20% level of 
the one from tttt scattering, which essentially refiects the experimentally observed K/tt ratio. Since tt-tt scattering 
could provide the dominant contribution to the overall low-energy yield (i.e., in the low-gt spectra of URHICs), the 
strangeness component may be significant. Also, the role of the kaon-chemical potential, ^k, is equally relevant for 
maintaining a substantial emission rate for n-K scattering toward lower temperatures (compare left and right panels 
of Fig. [6]) as is /x^r for tt-tt scattering. 

Additional sources of photon Bremsstrahlung in meson-meson scattering have been examined in Ref. [l2| within the 
soft photon approximation. At a temperature of T = 200 MeV and at low energies, qo < 0.5 GeV, Bremsstrahlung off 
TTTT scattering was found to exceed the one off ttK scattering by about a factor of ^5, roughly in line with our results 
(we find a factor of '^3-4). Furthermore, TTp —^ np^ processes were found to give a low-energy photon yield comparable 
to ttK —f ttKj, while other channels {KK, Kp, -kK*) are suppressed by an order of magnitude or more. However, 



8 

at lower temperatures, the rate for irp -^ 7rp7 will decrease faster than the one for ttK -^ ttKj, due to thermal 
suppression caused by the larger p-mass relative to the kaon mass (e.g., at T=120 MeV we estimate the relative 
suppression of former relative to the latter to be a factor of ^1.7, which includes the effects of the meson-chemical 
potentials; without the latter, the suppression will be larger), ttt] channels are suppressed relative to ttK, due to the 
mass and small degeneracy of the rj, which, in addition, is electrically neutral. Thus, for meson-meson Bremsstrahlung 
in heavy-ion collisions, one can expect all channels other than tttt and ttK to contribute on the ~20% level or less. 

Finally, let us compare the absolute magnitude of our rates to other calculations available in the literature. For 
T — 200 MeV, our tttt -^ ttttj rate for fj,^ = closely follows the result of Ref. |9i] (obtained in SPA), only that our 
rate is a factor 2 larger. This is precisely the enhancement we deduced from Fig. |4]as being due to the improvements 
in the present work. On the other hand, our results do not compare well with those found in Ref. [l2|. Fig. 3, where 
at go = 250 MeV and T = 200 MeV the rates are larger by at least a factor of ^2, despite being evaluated in SPA. 
The origin of this discrepancy deserves further study. 

IV. FURTHER UNCERTAINTIES IN THE EMISSION RATES 

Before we apply our Bremsstrahlung rate to estimate low-energy photon production in heavy-ion reactions at the 
SPS, we would like to illustrate two additional uncertainties in the determination of the emission rates. The first one 
relates to the well-known Landau-Pomeranchuk-Migdal (LPM) interferences, while the second one concerns possible 
medium effects on the meson exchanges which build up the meson-meson interactions. 

A. Landau-Pomeranchuk-Migdal Effect 

When the formation time of the emitted photon, Tform ^ l/^Oi becomes of the same order as the time between 
subsequent tt-tt collisions, r™", coherence effects are expected to become important (LPM effect [29|). For soft photon 
and dilepton production in a dense hadronic gas, the LPM interference has been studied in Refs. [3, [20,[2i|. Here, 
we put these estimates into context with our updated rate calculations. First, we recall the expression derived in 
Ref. ^(l| , where the LPM effect was implemented in SPA resulting in a differential cross section of photon production 
off TT-TT scattering as 

/da''\ 2a / 2 {1 - cos^O) \ ,^^. 

\d'^q / (27r)^ \ a^ + qQ[l — vcosO)'' / 

(•) indicates averaging over the time between successive collisions and over all the velocities after the collision (for 
simplicity, we choose an average velocity v = 0.5c), (T^tt denotes the elastic tt-tt cross section (which factorizes in SPA) 
and l/a = r™" ~ I/Ptt is the pion collision time (or inverse scattering width). The latter, in turn, depends on the 
hadron densities and pertinent pion cross sections, schematically given by 

r^ = ^nh{T,fih) a^h Vrei , (18) 

h 

where the summation is, in principle, over all hadrons in the heat bath (h — tt, K^ iV, A, . . .; VreX'- relative velocity be- 
tween TT and K). For a selfconsistently calculated pion gas at (T, /Itt) — (150, 0) MeV with realistic tttt interactions [30], 
one finds r™" ~ 6 fm/c, while a hot ttTVA gas representative for SPS conditions leads to t^°" ~ 2 fm/c [31[. In 
Fig. [7] we used both values to bracket the uncertainty in the suppression of low-energy photons, by plotting the ratio 
of Eq. P?)) to the incoherent limit (a = 0). For t^°" = 6 fm/c, the quenching factor is rather small reducing the 
photon emission rate by no more than 15% even at the smallest energies considered here, go = 100 MeV. This is 
consistent with Ref. [2l| where a slightly different coherence factor has been derived, amounting to a multiplication 
of the incoherent rate by [((7ot)^/(1 + (goT")^)]^; for t^°" ~ 6 fm/c and go = 100 MeV the latter turns out to be 0.81. 
However, for r™" = 2 fm/c and go = 100 MeV, collision and formation time are of comparable magnitude, and the 
quenching factor becomes substantial, between 0.4 and 0.25 according to the two ways of estimating it. 

B. Schematic Estimate of Medium Effects 

When hadrons are embedded into a hot and dense medium, their vacuum properties (mass and width) are expected 
to change, which is encoded in pertinent spectral functions, 

A.(fco, fc) = -2 ImA.(fc^) ^ -2 Im ( ,. _ ^. _^ ^^^(,^^ ,) ) , (19) 
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FIG. 7: Influence of the LPM efi'ect on thermal photon bremsstrahlung from tttt scattering at T = 150 MeV. 

where S/,(A;o, fc) is the self-energy of meson h as induced by interactions with the surrounding matter particles. For the 
present purpose, the relevant mesons are the p and the ^^a''' which figure into the meson-exchange propagators in the 
amplitudes (since we work in Born approximation, no explicit pion propagators appear). Dilepton spectra in heavy- ion 
collisions have revealed ample evidence that the in-medium p-meson spectral function is substantially modified. While 
the low-mass enhancement observed by CERES/NA45 in P6(158 AGeV)-Au collisions jS^] is compatible with both 
the assumption of a dropping p-mass and a broadening as predicted by hadronic many-body theory [23| , more recent 
NA60 dimuon spectra |33[ clearly favor the latter [34l| . Concerning the "cr" meson, 2-pion production experiments off 
nuclei |35l . |36[ indicate substantial medium effects as well, in terms of an enhancement near the two-pion threshold 
which increases with nuclear mass number A. Whether this is due to a dropping "cr" -meson mass [37[ or hadronic 
many-body effects [38l| (or both), is not clear at present. At finite temperature, the (the approach to) chiral symmetry 
restoration (implying degeneracy of tt- and u- spectral functions), is suggestive for a cr- mass approaching the pion 
mass [39J. It is therefore desirable to investigate to what extent Bremsstrahlung off tt-tt scattering is sensitive to 
modifications of "cr" and p mesons. Previous studies along these lines [40, |4l|, |42| have mostly addressed hadronic 
photon emission at higher energies (in connection with a dropping p mass). 

In light of the above, we here focus on the effects of a reduced cr-mass on low-energy photon emission. We will 
not address changes in the p-meson properties, as a simple dropping mass appears to be ruled out while a consistent 
implementation of changes in the width in the propagator requires accounting for the pertinent processes in the photon 
production rate as well (for the same reason, we do not consider width changes in the cr-propagator). 

When recalculating photon Bremsstrahlung off tt-tt scattering by reducing the mass of the a meson to m„ — 0.45 GeV 
(as compared to 0.525 GeV before), the thermal emission rate barely changes, cf. Fig.[H] This can be readily understood 
by recalling in the elastic tt-tt cross section the cr-meson contribution plays a minor role, cf. right panel in Fig. [2l 
Nevertheless, the smallness of the effect when dropping the a mass by ^^15% is not necessarily expected, but forces 
us to conclude that low-energy thermal photon spectra are not sensitive to (partial) chiral symmetry restoration in 
the scalar channel (for completeness we state that the rate of low-energy photon production increases by ^ 25% if 
one were to reduce the p-meson mass from its vacuum value of 0.77 GeV to 0.65 GeV). 

V. PHOTON SPECTRA IN HEAVY-ION COLLISIONS 

We now apply our Bremsstrahlung rates to the calculation of photon spectra at low transverse momentum in 
central P6(158 AGeV)-P6 collisions at SPS. We supplement the spectra with the baseline results of Ref. [11[^ and 



To avoid double counting, we have removed the contributions from p — > 7r7r7, tttt — > p7, K* — > Kix^ and txK — > K*^ from the rates 
of Ref. JJ|, as they are part of the tttt — > 7r7r7 and tttt — > 7r7r7 processes when using full tttt and vK cross sections; the latter include 
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FIG. 8; In-medium effects on thermal photon production implemented by a dropping "cr"-meson mass in Bremsstrahlung off 
TT-TT scattering. Dashed {ma = 0.525 GeV) and solid line {ma = 0.45 GeV) are almost indistinguishable due to the minor role 
of "cr" -exchange in the total elastic tt-tt cross section at low energy. 



compare the total to the recent WA98 data [15[. As in Ref. [Uj, the folding of the thermal emission rates over 
the space-time evolution is performed with the help of an expanding fireball evolution which has been parametrized 
to reflect basic features of hydrodynamic models. With a hadronic resonance-gas equation of state (with chemical 
freezeout at (/x^, T^h) = (250, 175) MeV) and a participant-nucleon number of Npart = 340 to represent the 10% most 
central collisions, the fireball volume at chemical freezeout is fixed to approximately reproduce the observed hadron 
multiplicities (after chemical freezeout, the evolution is constructed using entropy and baryon- number conservation, 
as well as meson-chemical potentials to preserve the correct hadron ratios). After averaging over the appropriate 
rapidity interval of the WA98 acceptance the thermal photon spectra follow as 



dN 



1 



dy / dTVFBir)qo 



dW 
~d?^ 



(20) 



with ijmin — 2.35 and Umax = 2.95 (Ay — 0.6). The fireball volume is taken to evolve with proper time as 

Vfb{t) 



-7r(zo + v^^^oT + 2«z'^^)(^T,o + -^arr'^f 



(21) 



with zq = 1.8 fm (corresponding to a formation time of 1 fm/c and translating into an average initial temperature 
of To = 205 MeV), rr,o = 6.25 fm, a^ = ar = 0.045 c^/fm, and w^ = 0.6c (note that only half of the total fireball 
volume must be used to reproduce the correct dNch/dy). QGP and mixed phase are completed after th — 4.46 fm/c 
and thermal freezeout is assumed at Tfo = 13 fm/c (plus one extra fm/c to allow for strong resonance decays). The 
numerically inferred temperature evolution in the hadronic phase can be conveniently parametrized as 



with ci=1.3, C2=0.639/fm, and a "radius" 



T{r) = 
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(unstable) p and K* mesons which should not be counted twice. 
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FIG. 9: Direct low-qt photon spectra as measured in central Pb-Pb collisions at SPS [la ] compared to thermal emission spectra 
from an expanding thermal fireball with QGP and hadronic phases [llj]. In the left panel (taken from Refs. [ly, [ij) QGP [431 ] 
(dash-dotted line) and hadron gas (HG) emission ^11] (dashed line, without meson-meson Bremsstrahlung) add up to the solid 
line, while the upper dashed line additionally includes an estimate [Ig] of Bremsstrahlung off tt-tt scattering in the soft-photon 
approximation (SPA) and without final-state Rose enhancement factors, representing the rate depicted by the solid line in 
Fig. 131 (the contributions from p —* -k-k^ decays in the HG rates of Ref. [ill ] have been removed to avoid double counting). The 
downward arrows represent 90% upper confidence limits set by WA98 in an earlier measurement [l4| . The curves in the right 
panel are as in the left panel except that the upper dashed line now includes the improved rates from elastic tt-tt and -n-K 
Bremsstrahlung (without SPA, including final-state Rose enhancement and with a similar caveat on double counting as in the 
left panel, cf. footnote (2} • 



Parameterizations of the numerically computed Bremsstrahlung rates along the temperature-chemical potential tra- 
jectory used here (as well as separately in goj T and /x^) are quoted in Appendix [Ul 

In Fig. [5| we summarize our results for direct photon spectra in central P6(158 AGeV)-P6 in comparison to the 
low-gt data of WA98 [l^ at the CERN-SPS. As noted in the introduction, these data are extracted via a photon 
HBT method requiring an assumption on the emission source size implying an extra (systematic) error of the data 
(in addition, strictly speaking, the weighted average of each data point does not correspond to the center of the 
transverse- momentum bin, but is shifted to somewhat smaller qt, due to the rather sharp increase of the yield toward 
lower qt). The left panel of Fig. [21 reproduces the previous calculations for QGP and hadron gas emission of Ref. [11| 
(lower three curves), which are in good agreement with earlier WA98 spectra at qt > 1.5 GeV, but significantly 
underestimate the low-momentum data at gt = 0.1 — 0.3 GeV. This is still true upon inclusion of Bremsstrahlung 
off TT-TT scattering [la . [l7| treated in SPA and without final-state Bose enhancement, even though this contribution 
is noticeable. However, if the improved Bremsstrahlung rates from tt-tt and w-K interactions are employed, the 
discrepancy is reduced appreciably, see right panel of Fig. [9] (note that the Bremsstrahlung yield decreases faster 
with energy than the other hadronic contributions, thereby not upsetting the agreement of the previous calculations 
with the experimental spectra (upper limits) at higher energies). These findings are presumably to be considered 
as an upper limit since at this point we do not account for LPM suppression effects as discussed in Sec. IIV Al The 
thermal low-gt photon yield is furthermore found to approximately scale with the fireball lifetime, which reiterates 
the importance of the later hadronic stages (on the other hand, the sensitivity to transverse flow effects is small). 

We also note that baryon-induced effects, which enter through a medium modified p spectral function [28] carried 
to the photon point and which dominate the contribution labeled "HG" at the depicted momenta, are appreciable (at 
the 30-60% level in the range qt = 0.1 — 0.3 GeV). Besides resonance decays (e.g., A,iV(1520) — > -^7), these include 
Bremsstrahlung from tt-N, N-N, and A^-A interactions (as following from ttNN~^ and ttAA^~^ excitations in the 
pion cloud of the p-meson), which have been constrained by (the inverse process of) photoabsorption on the nucleon 
and nuclei [44]. In Ref. 10], similar contributions were claimed to be smaller. 
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FIG. 10: Left panel: two-Gaussian fit (with centers at j/o = ±1.1 and width ay = 1.0) to the experimental proton rapidity 
distributions M7] in central P&(158 AGeV)-Pb. Right panel: coherent photon spectra from Pb-Pb collisions at SPS using final- 
state proton rapidity distributions from experiment (dashed line) or (5-functions centered at yp,T = ±0 (solid line), compared 
to the experimental data [l5i] . 

VI. COHERENT EMISSION FROM HEAVY-ION COLLISIONS 

Another source of soft direct photons in heavy-ion collisions that we would like to study here is (nonthermal) 
Bremsstrahlung emission due to coherent radiation off the protons within the incoming nuclei |45l . |46| . In the long- 
wave length limit, A^ — 2TT/q 3> D {q: photon 3-momentum, _D ~ 10 fm: typical system size), phase interferences are 
suppressed and the radiation can be assessed via the rapidity shift that the protons undergo from the initial nuclei 
to the finally observed distribution. Including both target and projectile, and for central collisions, one finds the 
following form of the photon spectrum [43, |4g| , 
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1 — UCOS0 1 — vpcos6 
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where vp is the beam velocity, y — tanh v is the final-state (FS) proton rapidity, and 



dN 
dy 



and 



dN 
dy 



(24) 



(25) 



are the pertinent FS rapidity distributions for projectile and target nuclei, respectively. The latter can be reasonably 
well fitted to experimental data [47| using two Gaussians, cf. left panel of Fig. [TOl In the right panel of Fig. [TUl we 
summarize our results for coherent emission from the longitudinal deceleration of the incoming protons. With realistic 
final rapidity distributions (dashed line) the photon yield constitutes only a few percent of the WA98 low-energy data, 
which presumably is an upper limit since already for go = 0.1 GeV, the pertinent photon wavelength, A-^ ~ 12 fm, is 
of the order of the system size rendering phase variations (interference) potentially relevant. Even with the extreme 
assumption of (5-function like proton rapidity distributions at yp,T — (full stopping), the coherent emission yield 
falls significantly short of the data. We thus conclude that this mechanism cannot be at the origin of the observed 
enhancement over previous predictions [ll|. 



VII. SUMMARY AND DISCUSSION 



In this article we have revisited the problem of soft photon production from hadronic matter which has recently 
received renewed interest by the observation of a large excess by the WA98 collaboration in central Pb-Pb collisions 
at the CERN-SPS. Earlier predictions using two-body mesonic and baryon-induced thermal photon production sig- 
nificantly underestimated these data. We have focused on a more elaborate re-evaluation of Bremsstrahlung from 
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meson-meson interactions. Based on a Ucm(l)-gauged meson-exchange model for elastic tt-tt and ir-K scattering (in 
Born approximation), our main improvement over previous calculations consists of a full numerical treatment of the 
collision integral which goes beyond the commonly employed soft-photon approximation and allows for the inclusion 
of final-state Bose enhancement factors. We have explicitly constructed pertinent contact interactions to preserve 
electromagnetic gauge invariance of the underlying Bremsstrahlung amplitude. The combined effect of the two im- 
provements amounts to an increase of the thermal tt-tt Bremsstrahlung rate by up to a factor of 2 at low photon 
energies (qq = 0.1 — 0.5 GeV), while tt-A' scattering adds another ^--^20% (essentially reflecting the K/tt ratio in the 
hadronic gas) . As a consequence, meson- meson Bremsstrahlung becomes the dominant thermal photon source at low 
energies. When combined with previous rate calculations, the convolution over a standard thermal fireball evolution 
reduces previously found discrepancies with the low-momentum WA98 data. This comparison, however, does not 
include the LPM effect, which we have estimated to suppress the soft emission rate by up to several tens of percent. 
We have furthermore shown that the Bremsstrahlung rates are rather insensitive to medium effects on the "cr" -meson 
in form of a reduced mass, due to the fact that the elastic tt-tt cross section is dominated by (s-channel) p exchange. 
We finally studied the possibility of coherent Bremsstrahlung induced by rapidity shifts of the incoming projectile 
and target protons. When using realistic final-state distributions, we have found this contribution to be essentially 
negligible. 

While our investigations suggest that thermal emission from hadronic matter might be at the origin of the low-qt 
WA98 enhancement, more work is required to quantify the theoretical estimates, e.g. by merging medium effects on the 
p-meson (as implicit in rate estimates using in-medium spectral functions) with the newly computed Bremsstrahlung 
rates, as well as by explicit implementations of the LPM effect. In addition, an experimental confirmation of the 
observed enhancement, e.g. at the Relativistic Heavy- Ion Collider (RHIC), would be very illuminating. 
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APPENDIX A: EVALUATION OF THE BREMSSTRAHLUNG PHASE SPACE INTEGRAL 



In this appendix, we derive the expressions for the photon emission rate as used in Sec. IIIII Starting from the 
kinetic theory expression, Eq. ([7]) in the main text, for tt + tt^tt + tt + 7, 

dR-* ^^ f d^pa d^pt, d^pi d^P2 



d^q J 2Ea{2T:f 2Ei,(2T:f 2Ei{2TTf 2E2{2ttY 

•k{2tt) 6 [pa+Pb-pi-p2-q)\Mi\ 2{2tt)^ ' ^ ' 

we perform the integrations over the energy-momentum conserving delta function to obtain 
dW N f d^pa d^pb f d^pi 



qo 



■/^f^/fiAMi + /3][i + /4 



d^ 16(2TT)iiy Ea Eb 

X \M^\'^6{{Ea + Eb-Ei- qaf - (pa + Ph - Pi - q)^ - ml) 

^^ J \Pa\dEadna J \pb\dE,dn, J \p,\dE,dnjj2[l + /3][1 + /J 



16(2tt)i 

X \M^\^5{{Ea + Eb-Ei- qof - {pa + Ph - Pi - q)' - m^) . (A2) 

We evaluate this equation by the Monte Carlo method, choosing the coordinates for each particle as follows, 

Pa = (-Ba,|Pa|sin6'aCOS0a, |Pa|sin6'aSin<?!)a,|Pa|cOs6'a), 

Pb = (-B6,|Pb|sin6'bCOS0b, |Pb|sin6'bSin0b,|Pb|cos6'b), 

Pi = (-Bi,|pi|sin6'icos0i,|pi|sin6'isin(/)i,|pi|cos6'i), 

q = (<zo,0,0,|q|), 

P2 = Pa + Pb - Pi - q , (A3) 
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with 

Thus, Eq. (jA2|) can be rewritten as 



E^ = Xr^i+Pi i^ = a-,b,l,2,q) 



(A4) 



«o 



d^q 16(27r) 



^ f \Pa\dEaSine,d0,d(f>a f \pb\dEbSine^debd(f>b f \pi\dEism0ideid(f,i hf 2[l + /3][1 + h] 



x\Mif5{{Ea + Eb-Ei -qof - (Pa+Pfc-Pi -qf -ml 



(A5) 



In the argument of the delta function, one has 

(Pa + Pb - Pi - q)^ = |PaP + 2Pq • pb - 2pa • Pi - 2pa • q + (pfc - Pi - q)^ 

= 2|pa||ph|sin6'aSin6'bC0S(/)aC0S(/)b + 2|pa||pfc|sin0aSin6'bSin(/)aSin(/)b + 2|pa||p6|cos6'aCOs6'b 
-2|pQ||pi|sin6'aSin6'iCOS(/)aCOS0i - 2|pa||pi|sin6laSin6lisin0aSin(?!)i - 2|pa||pi|cos6laCOs6'i 
-2|p,||q|cos0a + |Pa|' + {Pb - Pi - q)' . (A6) 

Upon integrating over 9a, we arrive at 

^^7 AT r r r Ia.«.|2 



90 



d^q 



— ^^ J \pa\dEadcf>a J \pb\dEbSm9^d9bd4>^ J \pi\dEism9id9id4>ifif2[l + /3][1 + M^^ (A7) 



with 



(A8) 



^=|^'(cos0DI- 
In the above, cos^J^ is the root of function i^(cos0a), 

V3(cos6'q) = (2|po||pb|sin6'bCos0aCOS(/)b + 2|pa||Pb|sin6'bSin0aSin0b - 2|pa||pi|sin6'iCOS(/)aCOS0i 
-2|pa I |pi |sin6'isin(/)asin</)i ) yT~~cos2^ 
+ (2|pa||pb|cos6'b - 2|pq||pi|cos6Ii - 2|pQ||q|)cos0a 
+ |Pa|' + (Ph - Pi - q)' - [Ea + Eb-E,- q^f + m2 = 0. (A9) 

with 

(Pb - Pi - q)^ = |Ph|^ + |PiP + Ip-jI^ - 2|p6||pi|sin6'bsin6'icos0bcos(/>i - 2|p6||pi|sin6'bsin6'isin</)bsin</ii 

-2|pb||pi|cos6'bCos6'i - 2|pb||q|cos6'b + IpiUqlcos^i. (AlO) 

There are two additional constraints on the eight integration variables in Eq. (IA7|) : 

a. For any configuration of the eight variables, Eq. (jA9p must have at least one root. 

b. Ea + Eb > El + qa + m2. 



APPENDIX B: GAUGE INVARIANCE OF PHOTON PRODUCTION AMPLITUDES 

In this appendix, we demonstrate gauge invariance for the process tt+tt^ — > 7r+7r^7. In this case the charges are, 
Qa — Qi = 1, Qb — Q2 ~ —1, and there is no w-channel contribution. To verify A^'' ■ q^ ~ for the i-channel 
amplitude, we have 

M'^ ■ q^, = e{J^Mt{pa-q,Pb,Pi,P2) + JbMt{pa,Pb-q,Pi,P2) 

+ j'^Mt{Pa,Pb,Pl + q,P2) + J2-Mt{Pa,Pb,Pl,P2 + <?)) ' Q,, 
= e{-Mt{Pa ~ q,Pb,Pl,P2) + Mt{Pa,Pb - Pa,Pl,P2) + Mt{Pa,Pb,Pl +q,P2) " Mt{Pa, Pb, Pi, P2 + q) 

{Pa-Pl- q- Pl)Pb ■ P2 {Pa -Pl+Pa- q)Pb ' P2 



Aegl^^^F\pb-p2f] 

+ F[{pa--piY] 



[Pb - P2)^ ~ml + irriaVa [pb - P2f - ml + im^T 

Pa ■ Pl{Pb ■ P2 - q ■ P2) , Pa ■ Pl{Pb ■ P2 + Pb ■ q) 



[Pa -PiY -ml + irn^T^ {pa - PiY - ml + im^T^ 



Wa 



p2r^_ „ n2i -iq-Pl+q-Pa)Pb-P2 _^_p2u _ n2i Pa ' Pl{q ' P2 + q ' Pb) 



F'[{pb~P2) 



[pb - P2Y - ml + imaVa 



{jpa ~ piY - ml + imaVo 



(Bl) 
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so that the pertinent contact term is deduced as 



M['^'-q, = F^[ip,-p,r] 



4e5^^^ {Pa-q+Pi- q)Pb ■ V2 ^2 



{Pb - P2Y - ml + im^Ta 



F'[{pa-Pi) 



1 4egg^^Pa ■ PijPb ■ q + P2 ■ q) 
[Pa -piY -ml + im^Ta- 



resulting in A^^ • g^ + A^^ ' • g^ = 0. Similarly, for the s-channel, we have 



4eff; 



"''^ [{Pi +^2)^ - ml + imaTa {pi + P2Y - ml + imcTa 

Pa ■ Pbipi ■ P2 + q ■ P2) , Pa ■ Pb{Pl ■ P2 + Pi ■ q) 



{Pa + PbY - ^2 + irriaTa (Pa + PbY " ml + im^Ta 



Wa 



-{q-Pb-q-Pa)pi ■P2 



Pa -Pbiq-Pi - q-P2) 



_ {Pi +P2)'^ -ml + im„Ta {pa + PbY - ml + im^Ta 
and the required contact term follows as 

■^(^al-KTriPb -q-Pa- q)Pl ■ P2 Aegl^^Pa ■ Pb{Pl ' q - P2 ' ?) 



TVjWp 



9m 



{Pi + P2Y -ml + im^Ta {pa + PbY - ml + im„T„ 



<(c)m 



(B2) 



Mi^-qf, = e{J^MsiPa ~ q,Pb,Pi,P2) + JbMsiPa.Pb- q,Pi,P2) 

+ J^Ms{Pa,Pb,Pl + q,P2) + J2Ms{Pa,Pb,Pl,P2 + ?)) ' 9^ 
= e{~Ms{Pa - q,Pb,Pl,P2) +Ms{Pa,Pb -Pa,Pl,P2) +Ms[Pa,Pb,Pl +q,P2) ~ Ms[Pa,Pb,Pl,P2 + q) 
[Pa-Pb- q- Pb)Pl ■ P2 {Pa -Pb-Pa- q)Pl ' P2 



(B3) 



(B4) 



Again, we arrive at M.'^ ■ q^ + AAl ■ qfj. — 0. It is straightforward to test the rho-meson exchange channels along 
similar lines. 



APPENDIX C: RATE PARAMETERIZATIONS 



The thermal production rate of photons from Bremsstrahlung off tt-tt and ir-K scattering as employed in Sec. |V] 
can be conveniently parametrized as 



dR^ 



-0.00026421 



90- 



0.0075394 
0. 13007 Y^Tiyi 



exp 



-2(T-0. 28351) 
0.13007^ 



<Pq 



-0.00018184- 



0.01277 

0.00073983 

7r[4(go -0.054587)^+0.099656^1 



90- 



dR^ 



[1. + 104(0.18 -T)i-4(qo- 0.1)] 
-0.0000411907 -L 0.0013114 



K 



0.12037-v/7r/2 



exp 



2(T-0.27292)"' 
0.120372 



d3g 



-0.000021454 



0.00227 

0.0001459 

7r[4(go-0. 043159)2+0. 0856682 



[l.-f 605(0.18 -r)i-7(go- 0.1)] 



(CI) 



These expressions reproduce the exact rates within ~3% in the energy interval between 0.1 GeV and 0.5 GeV and in 
the temperature range 0.1-0.18 GeV. Pion and kaon chemical potentials are included as ii-k{T) and ij.k{T) according 
to the fireball evolution employed in this work for central P&(158 AGeV)-P6 collisions, cf. eq. (|16p . Temperature T 
and energy go are in units of GeV, and the resulting rates are in units of l/GeV^fm'*. For more general use, we also 
provide parameterizations which explicitly resolve the qo, T and fi^^ dependencies: 



90 



dR^^ _ [l-Fgo/(0.011-F0.56T-0.565T2)]-4(o.92-|-0.17e^"/°°247) 



90- 



d^q 

dR-nK 
1^ 



3(0.87 + 519e-^/oo26)(o.92 -f 372e-^/0-026) 



[1 + go/(-0.024 + 0.773T)]-4(o.8 + 0.223e^'^/°°53i4yo.92 + 0.176^^-/0 
(14.87 -I- 3856e-'r/"032)^_92 + 372e-'rAi^026(0.89 -f 143e-^/oo26) 



0247 



(C2) 
(C3) 
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These expressions reproduce the exact rates within ~5% in the energy interval between 0.1 GeV and 0.5 GeV and in 
the temperature range 0.1-0.18 GeV. 
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